A Model for Conductive Percolation in Ordered Nanowire Arrays 
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The combined processes of anodization and electrodeposition lead to highly ordered arrays of cylin- 
drical nanowires. This template-based self-assembly fabrication method yields nanowires embedded 
in alumina. Commonly, chemical etching is used to remove the alumina and free the nanowires. 
However, it has been experimentally observed during the etching process that the nanowires tend 
to form clumps. In this work, the nanowires are modeled as elastic rods subject to surface interac- 
tion forces. The dynamics of the model give rise to the aforementioned clumping behavior which 
is studied via percolation theory. This work finds that percolation takes place with probability 
P ^ {t — tc)^ , where the exponent x — 2.8 and tc is the time at which percolation takes place. 
The critical exponents which entirely determine the system are found to be for (dimension) d — 2, 
/3 = 2.1, 7 = 0.57, A = 2.7, a = -2.8, u = 2.4, and 6 = 1.3. 

PACS numbers: 64.60.ah, 61.46.-w, 02.60.Cb 



I. INTRODUCTION 

The fabrication of nanowires through anodization and 
electrodeposition provides an interesting opportunity for 
the study of percolation. Before presenting the details of 
how this comes about, we first introduce the related key 
topics. 

A. Template Based Fabrication 

Under proper conditions, the anodization of aluminum 
creates a thin layer of alumina containing a highly or- 
dered hep array of pores [1]. These pores, though typi- 
cally of a uniform size in a well made sample, can have 
diameters ranging from 10 nm to 100 nm. Their length, 
being a function of time, can easily extend to the micron 
scale. Electrodeposition of metals such as Au, Ag, Co, 
Ni, or Fe results in a self-assembled and highly ordered 
array of cylindrical nanowires. 

Post-fabrication, the nanowires exist as a bulk sample; 
they remain embedded within the alumina. If it is so 
desired, a chemical etch can remove the alumina and re- 
lease the nanowires into solution. Figure [T] illustrates the 
etching of a sample at different times. 

B. Conductive Percolation 

Percolation theory has been used to study conductivity 
in a variety of electrical systems j2|, |3| , but more recently 
has seen applications to carbon nanotubes 0, ^ and the 
growth of V2O5 nanowires 0. In these examples, the 
nanostructures are initially suspended in solution and al- 
lowed to settle on the surface of a substrate. The re- 
sulting random network is probed for a conductive path 
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FIG. 1: Seen from a cross sectional slice, the anodized alu- 
mina acts as a template for nanowire growth. The sequence 
A-D illustrates various etch times where A is no etching; B is 
moderate etching which exposes just the tips; C is extended 
etching which results in clumps of nanowires; and D is total 
etching which releases the nanowires. Region C is of partic- 
ular interest because this is when the percolation phenomena 
begins to takes place. 

as a function of surface coverage. Generally, it is found 
that at some critical density, the conductivity rapidly in- 
creases as more electrical paths become available for the 
current. 



C. Clumping of Nanowires 

The key observation which relates percolation and 
nanowire fabrication in alumina templates, is the behav- 
ior of clumping. As has been reported III, [8|, during 
the etching process, nanowires tend to clump into bun- 
dles. If the process is interrupted at various times, the 
average clump size (i.e. number of nanowires in direct 
physical contact) is seen to grow. Figure [2] shows a typi- 
cal example of this behavior. 

In an experimentally realizable scenario, the clumped 
nanowires can be tested for conductive percolation by 
placing electrical contacts along two opposing edges of 
the sample. When the template has been sufficiently 
etched (as in figure [T]G), an electrical path should be ob- 
served in the form of a percolating cluster. 

Inspired by these results, we consider a two dimen- 
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FIG. 3: In calculating the van der Waals interaction energy, 
we use the geometry of two cylinders of equal length L, sepa- 
ration radius R\ and R2. The vectors in Eqs.JS]) are defined 
as shown. 



FIG. 2: During extended etch times, the exposed portions of 
the nanowires bend from their initial position and clump into 
groups. This results in large clusters of nanowires making 
physical contact and leads to the possibility of conductive 
percolation. In this SEM image, the etching has taken place 
for sufficiently long time. Clumps have begun to take shape 
as the nanowires bend due to surface forces. 



sional array of nanowires, modeled as elastic rods, and 
subject to van der Waals interaction forces. In the model, 
we consider the nanowire length a linear function of time 
to simulate the effect of etching. The dynamics of this 
system are simulated in MATLAB and it is observed that 
given sufficient etch time, the nanowires bend in the di- 
rection of their neighbors. This gives rise to clusters of 
variable size which we interpret from the view of perco- 
lation theory. In section [III the van der Waals interac- 
tion force is derived for the cylindrical geometry of two 
nanowires. In section Hill we review the equations of equi- 
librium which describe the bending of elastic rods. In 
section |IVl the numerical results are presented and dis- 
cussed in terms of percolation. In section IVl two experi- 
ments are proposed which would quantitatively verify the 
model described in this work. 



II. SURFACE INTERACTIONS 

To calculate the interaction forces between two macro- 
scopic bodies, we use the non-retarded van der Waals 
potential. If we consider the bodies to be composed of 
discrete point-like atoms, the total energy of interaction 
can be written as the sum 
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We denote by r^j the separation between two given par- 
ticles, the i^^ in the first body and the j*^ in the second. 
The constant C is the coefficient of the atom-atom pair 



potential. Passing to the continuum limit. 
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where g is the number of atoms per unit volume. Because 
we are interested in the interaction of two nanowires, the 
natural coordinates are cylindrical (see figure [3]). It can 
be seen from simple vector algebra that 

D = (i?i+i?2+c^,0,0), 

ri = (picos6>i,pisinl9i,2;i) , 

1*2 = {p2 COS O2 , p2 sin 6>2 , 2^2 ) , and 

r = D + r2-ri. (3) 

Expanding the denominator of Eq.([2|) in terms of com- 
ponents, the integral can be calculated. The final result 

is m 
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where A is the Hamaker constant defined by 
A = ti'^CqiQ2- 
The force arising from this potential is, 

^ ^ ""SV^ |d|5/2 \Ri^R2 
where n = d/d. Letting Ri = R2 = R^ 
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When immersed in aqueous solution, the Hamaker con- 
stant is not accurately given by Eq.([5|). Because interac- 
tions in the medium effect the pair-wise potential, exper- 
imentally determined results for metals in water must be 
used [10]. 



III. EQUILIBRIUM OF A BENDING ROD 

We now review the equations of equilibrium for a rod 
subject to external forces. The discussion closely follows 
pj| and only the main results will be presented. 

Consider a thin rod of circular cross-section. To an 
infinitesimal length, assign local coordinates ^, 7^, ( so the 
("-axis is aligned parallel to the rod. Furthermore, let r = 
r(/) be a curve parameterized by arc length describing the 
local coordinates, and hence, the unit tangent vector is 
given by t = dr/dl. We define 

n = d^/di, (8) 

to be the rotation of the local coordinates along the rod 
and notice that in the absence of torsion, t • 11 = r^<^ = 0. 
Moving along an infinitesimal length, the change in the 
tangent vector is given by 

dt d^ ^ 

Taking the cross product of both sides with t and rear- 
ranging, we find the expression 

n = tx^^t{t-Q). (10) 

In the case of pure bending, the last term is zero due 
to the absence of torsion. Let the moment vector M be 
defined by 

m = Em, (11) 

where E is Young's modulus and / is the moment of 
inertia of the cross section for a cylinder. As a condition 
for equilibrium, the total moment of the forces is zero, 
i.e., dlsA + (ir X F = 0. Dividing by dl, 

d'M/d\ = Y xt. (12) 

Rearranging Eqs.([TOl- [T2|), we find the equation of equi- 
librium for pure bending in a circular rod: 

^ dr ^^dr d^r , 

IV. NUMERICAL RESULTS 

To estimate the magnitude of the force in Eq.([7j) and 
the amount of displacement the tip of a nanowire will 
undergo, characteristic values are used p]. Isl. fiol. fl2l. . 
For Au nanowires with radius R = 20 nm, interpore sep- 
aration d = 10 nm, length L = 500 nm, and Hamaker 
constant A = 10~^^ J, the force F = 44 pN. The amount 
of defiection can be found by treating the nanowire as 
a cantilever with deflection Ad = FL^/3EI. For Au, 
Young's modulus E = 7.8 x 10^ Pa, the moment of iner- 
tia of the cross-section is / = {7r/4)R^ = 1.3 x 10~^^ m"^, 
and consequently, the deflection is Ad —19 nm. This 
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FIG. 4: A 5 x 5 grid is simulated for tmax = 3 to illustrate the 
results of the model. As shown, some nanowires are already 
bending to make contact with their neighbors while others 
remain unaffected. For computational reasons, the perimeter 
is considered infinitely rigid and those wires do not bend. 



estimate indicates that van der Waals forces are suffi- 
ciently capable of bending the nanowires. Consequently, 
we proceed to a more sophisticated numerical solution of 
the equations involving an array of nanowires. 

The array of nanowires is simulated using a 4*^ or- 
der Runge-Kutta method to solve Eq.([T3j) subject to 
the forces given by Eq.([7|) on a square grid. Each grid 
site {ij} corresponds to an individual nanowire and the 
forces acting on it are determined by nearest neighbor 
interactions solely concentrated at the tip. Because each 
nanowire is treated as an elastic rod, the net forces cause 
bending accordingly. The continuous etching of alumina 
shown in ffgure [T] is simulated by making the length of 
each nanowire a linear function of time. Thus, taking 
Zij = Vot where t G (0, tmax 

] and is constant, the pro- 
cess can be stopped at each unit of time and checked for 
percolation. Furthermore, the length of each nanowire 
is randomized to within 15% of the maximum value to 
reffect experimental observations. The boundries were 
treated as inffnitely stiff rods and hence, not subject to 
interaction forces. Figure [4] shows a 5 x 5 sample simu- 
lation where some nanowires have already made physical 
contact with their neighbors. 

In a 2D representation, figure [5] shows typical results 
for repeated simulations of various etch times. The blue 
line segments represent physical contact between two grid 
points and the largest clusters have been colored red. 
Each simulation was deemed percolating when a cluster 
spanned the full horizontal distance. Running 500 sim- 
ulations for each tmax = 1,2, ...,13, the probability of 
percolation was determined. Figure [6] displays these re- 
sults along with the fitted curve P (x {t — tcY ^ where 

= 4 and the critical exponent x = 2.8. This nonlinear 
scaling is indicative of the percolation threshold. 

Near the threshold, percolation theory predicts that 
the probability of percolation as a function of the prob- 
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FIG. 5: Typical results for simulations of various times. (A) 
t = 3; (B) t = 6; (C) t = 8; and (D) t = 13. Blue line 
segments indicate two nanowires that are touching. Blank 
sites indicate free-standing nanowires, and the largest clus- 
ter in each figure was colored red to illustrate the growth in 
cluster size as a function of time. Both (C) and (D) contain 
percolating clusters. 
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FIG. 6: The probability of percolation is plotted against the 
simulation etch time. It is seen that at t ~ 4 the probability 
grows like P oc (t—tcY^ where tc = 4 and the critical exponent 
X 2.8. 
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FIG. 7: The probability of percolation is plotted against the 
probability of physical contact between two wires. This scales 
like (p — Pc)^,P ^ , which aids in determining the full set 
of critical exponents. From the fit, it is observed that /3 = 2.1. 
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FIG. 8: The cluster size is plotted against the probability of 
physical contact between two wires. The singularity in the 
neighborhood of p ~ 0.4, is due to the presence of large clus- 
ters. As the threshold is crossed and the clusters become infi- 
nite in size, they make no contributions to the average cluster 



size and consequently, the curve decays as {pc —p) 
From the fit, it is observed that 7 = 0.57. 



determined by the equations 



ability of physical contact, P(p), will obey 

P(x(p-Pc)^, P 



(14) 



where jS can be related to other critical exponents of the 
system. Similarly, the mean cluster size as a function of 
the probability of physical contact, S{p)^ goes like 



5' (X {Pc-pY 



P^Pc 



(15) 



Prom figures [7| and IHl /3 = 2.1 and 7 = 0.57, can be read 
from the curve fits. The fuU set of critical exponents is 
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-13, 
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Therefore, in (dimension) d = 2 for the percolating 
nanowire network subject to van der Waals interac- 
tion, the complete set of critical exponents is given by 
(3 = 2.1, 7 = 0.57, and with the aid of Eqs.([T6]), 
A = 2.7, a -2.8, u 2.4, and S 1.3. 
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V. CONCLUSION 

The goal of this work was to quantitatively explain 
and describe the phenomenon of nanowire clumping dur- 
ing the etching of ordered nanowire arrays. Having ac- 
complished this task, we now outline two potential ex- 
periments which could be used to verify the predicted 
critical exponents. In the first case, a number of samples 
could be fabricated in the laboratory and by means of 
counting small subregions^ all the numerically simulated 
results could be readily verified. 

The alternative experiment involves measuring the 
conductivity across an array of nanowires while the alu- 
mina is etched. At time t = 0, all the nanowires would 
be embedded in the alumina (an insulator) and no cur- 
rent would fiow. At some time t = tc^ the first perco- 
lating cluster would form and the current density would 
increase from 0. Continuing the measurements for longer 
time, the conductivity would increase further as more 
pathways through the nanowire array opened up allow- 



ing a greater current to fiow. As t ^ oo, the conductiv- 
ity would approach a constant value when the clusters 
have fully formed and no more potential pathways exist. 
By means of the usual analysis in conductive percolation 
[14], the relevant statistical distributions can be deter- 
mined from this data. 

In conclusion, the qualitative observation of clumping 
was modeled by elastic rods subject to van der Waals in- 
teraction forces. Simulations in MATLAB revealed a per- 
colative phase transition, from which, the critical expo- 
nents were extracted. They are /3 = 2.1, 7 = 0.57, A = 
2.7, a = -2.8, ly = 2.4, and S = 1.3. 
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